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Orbits in a three-dimensional potential subjected to periodic driving, V{x'^,t) = (1 + 
mt,smLjt)Vo(x^), divide naturally into two types, regular and chaotic, between which transitions 
are seemingly impossible. The chaotic orbits divide in turn into two types, apparently separated 
by entropy barriers, namely 'sticky' orbit segments, which are 'locked' to the driving frequency and 
exhibit little systematic energy diffusion, and 'wildly' chaotic segments, which are not so locked and 
can exhibit significant energy diffusion. Attention focuses on how the relative abundance of these 
different orbit types and the transition rate between sticky and wildly chaotic orbits depends on 
amplitude mo, and on the extent to which these quantities can be altered by weak friction and/or 
noise and by pseudo-random variations in the driving frequency, idealized as an Ornstein-Uhlenbeck 
process with (in general) nonzero autocorrelation time tc- When, in the absence of perturbations, 
there exist large measures of both regular and chaotic orbits, the primary effect of weak noise is 
to increase the relative measure of chaotic orbits. Alternatively, when almost all the orbits are 
already chaotic, noise serves primarily to accelerate transitions from sticky to wildly chaotic behav- 
ior. The presence or absence of friction is unimportant and the details of the noise seem largely 
immaterial. In particular, there is only a weak, roughly logarithmic dependence on amplitude, and 
additive and multiplicative noise typically have virtually identical effects. By contrast, allowing for 
high frequency 'noisy' variations in ui tends to weaken effects of the driving, decreasing the relative 
measure of chaotic orbits and suppressing large scale energy diffusion. For both ordinary noise and 
noisy frequencies, the largest effects are observed for tc — » 0, the efficacy of the noise decreasing 
with increasing tc- 
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I. INTRODUCTION AND OVERVIEW 

Rapid energy diffusion in nearly periodic, oscillating 
potentials is an important phenomenon arising in a va- 
riety of physical settings, including collective relaxation 
in many-body systems interacting via long-range forces. 
One example is 'violent relaxation' P, the process 
whereby a galaxy involved in a collision or close en- 
counter with another galaxy evolves towards a (quasi-) 
equilibrium. Another is halo formation in an accelera- 
tor beam j^, where particles in an initially concentrated 
bunch are ejected to relatively large distances from the 
center. 

This efficient energy diffusion can often be attributed 
to parametric resonances between the frequency of the 
periodic driving and the natural frequencies of the or- 
bits {e.g., 01) which, for appropriate driving frequen- 
cies, tend to make many otherwise regular orbits strongly 
chaotic 1^. However, not all chaotic orbits exhibit 
systematic energy diffusion. In at least some time- 
dependent spherically symmetric systems chaotic orbit 
segments divide empirically into two types 0. 
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On the one hand, there are 'sticky' chaotic segments 
which are 'locked' to harmonics of the driving frequency 
and hence, albeit chaotic with positive finite time Lya- 
punov exponents, do not exhibit systematic drifts in en- 
ergy. On the other, there are 'wildly' chaotic segments 
which are not so locked and, as such, can and do drift in 
such a fashion as to exhibit large-scale energy diffusion. 
However, this distinction between 'sticky' and 'wildly' 
chaotic is not absolute. If, e.g., a sticky orbit is inte- 
grated long enough, it can become unstuck and start be- 
having as a wildly chaotic orbit. Nevertheless, the time 
required for such a transition can be quite long, hundreds 
of orbital times Id or more. 

This transitional behavior appears analogous, at least 
qualitatively, to diffusion through cantori in two-degree- 
of-freedom Hamiltonian system. Cantori act as 'entropy 
barriers' which impede, but do not prevent, phase space 
transport However, barrier penetration in such sys- 
tems can be dramatically accelerated by allowing for 
noise sj'o'l which, by wiggling the orbits, helps them to 
'find' suitable holes. 

In this setting noise too acts via a resonant coupling 
. White noise, which has zero autocorrelation time 
and power at all frequencies, tends to be comparatively 
efficient as a source of accelerated phase space transport. 
By contrast, colored noise, with a finite autocorrelation 
time, tends to be efficient only if the autocorrelation time 
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is less than or comparable to the orbital time scale, so 
that the noise has substantial power at frequencies for 
which the orbits have power. 

An obvious question, therefore, is whether transi- 
tions between sticky and wildly chaotic orbits in time- 
dependent Hamiltonian systems can also be accelerated 
by noisy perturbations. Here there are (at least) two dif- 
ferent types of 'noise' to consider. There is of course the 
possibility of 'ordinary' noise, either intrinsic or extrinsic 
(see, e.5.:[lll), incorporated simply as an extra term in 
the equations of motion, just as for a time-independent 
Hamiltonian system. However, there is also the possi- 
bility of allowing for a 'noisy' driving frequency, i.e., al- 
lowing the frequency to vary in a (near-)random fashion. 
Even a modest 'wiggling' of the frequency might weaken 
the resonant couplings, thus decreasing the effects of pe- 
riodic driving. 

In all this, one is interested primarily in the effects 
of low-amplitude perturbations. After all, it is obvi- 
ous that, for sufficiently large perturbations, solutions 
to the perturbed problem will be very different from so- 
lutions to the unperturbed problem. The real question is 
whether relatively small perturbations, which act in 'real ' 
systems but which one might naively exclude in physical 
modeling, are sufficiently important that they cannot be 
safely ignored. For example, even relatively low ampli- 
tude oscillations in the bulk potential associated with a 
charge-particle beam can prove important by facilitating 
the ejection of charges to an extended outer halo; and al- 
lowing for modest frequency variations can actually make 
the effect more pronounced ^bi- 
section H focuses on the empirical distinctions between 
sticky and wildly chaotic orbits in a strictly periodic po- 
tential, addressing, in particular, the issue of robustness. 
Earlier work established that these distinctions exist in 
spherically symmetric potentials. However, the situation 
is less clear for more complex potentials. What happens 
if, for instance, chaotic orbits exist even in the absence of 
the periodic driving? Could this tend to decrease the rel- 
ative measure of sticky orbit segments or, perhaps, elim- 
inate them altogether? Section HI focuses on the role of 
'ordinary' noise, modeled as an Ornstein-Uhlenbeck pro- 
cess; and Section IV considers the role of random vari- 
ations in the driving frequency. Section V summarizes 
the principal conclusions and then comments on poten- 
tial implications. 



II. CHAOS AND ENERGY DIFFUSION 

The objective here is to show that absolute distinctions 
between regular and chaotic orbits and short-time dis- 
tinctions between 'sticky' and 'wildly' chaotic orbit seg- 
ments, observed @ in time-dependent one-dimensional, 
6.17., spherically symmetric, potentials persist in more re- 
alistic time-dependent three-dimensional potentials; and 
to investigate for several examples how the relative mea- 
sures of different orbit types depend on properties of the 



potential. 

Attention focuses on time-dependences of the form 

y(x\<) = (1 + mosintJoi) X ^0(2:'), (1) 

with Vq an unperturbed time-independent potential. 
Two classes of potentials Vb were considered, namely 

^0(^,0 = -(Y^^ (2) 

and 

Vo[s,t) = --^, (3) 
1 + s 

each with 

s^ = {x/af + {y/bf + {z/cf. (4) 

For spherical systems, with a = b = c = 1, these reduce 
to integrable potentials well known from galactic astron- 
omy. The former corresponds via the Poisson equation to 
a so-called Plummer density distribution with a smooth 
central core (see, e.g., 0); the latter yields a 7 = 1 
cuspy Dehnen density distribution similar to that 
observed in many galaxies [l^ . The nonspherical gen- 
eralizations of the Plummer potential tend not to admit 
many chaotic orbits but, for appropriate choices of axis 
ratios, the nonspherical Dehnen potentials do admit sig- 
nificant measures of chaotic orbits. For the intermediate 
energies E —0.8 used for the numerical experiments de- 
scribed in this paper, a typical orbital time scale tjj ~ 5. 

Ensembles of > 800 initial conditions were generated 
by uniformly sampling an (unperturbed) constant en- 
ergy hypersurface, and these initial conditions were then 
evolved for a time t > 4096, allowing for a driving fre- 
quency LOfj so chosen as to trigger parametric resonance. 
The integrations also solved simultaneously for the evo- 
lution of a linearized perturbation, thus facilitating an 
estimate of the largest finite time Lyapunov exponent 
X. Orbital data recorded at fixed intervals were Fourier- 
transformed and the resulting data analyzed to establish 
correlations between the value of the largest Lyapunov 
exponent and the frequencies for which the x-, y-, and 
z-coordinates of the orbits have the most power. 

For all the models that were tested, involving pulsa- 
tions of both integrable and nonintegrable potentials Vq, 
the orbits divide naturally into two distinct types, namely 
regular and chaotic, between which transitions appear 
impossible. Estimates of the largest Lyapunov exponent 
for the regular orbits decay towards zero in the usual way 
(see, e.g.. U^): for chaotic orbits they remain strictly 
positive |l5j| . Moreover, regular orbits appear multiply 
periodic, whereas the Fourier spectra of chaotic orbits 
are broader band. 

It also appears that chaotic orbit segments subdivide 
into 'sticky' and 'wildly' chaotic. Sticky chaotic segments 
are clearly chaotic in the sense that they have positive fi- 
nite time Lyapunov exponents and more complex Fourier 
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spectra than do regular orbits. However, the spectra do 
exhibit striking regularities. In particular, the frequen- 
cies for which the transforms of the spatial coordinates 
have the most power are almost exactly equal to a har- 
monic of the driving frequency luq . In many cases almost 
all the power is at that special frequency; in others there 
is substantial power at nearby, mostly lower, frequencies 
but the driving frequency is still clearly dominant. That 
the orbits appear 'locked' to a harmonic might suggest 
that their energies cannot change significantly, a predic- 
tion which is readily confirmed. Plots of E{t) for sticky 
orbit segments often look nearly periodic visually, and 
the associated Fourier spectra tend to be very 

sharply peaked (although there is always some 'scruff' in- 
dicating that the energy is not exactly periodic). Wildly 
chaotic orbits have substantially more complex spectra, 
and they are not obviously locked to any harmonic of 
the driving frequency. Rather, as viewed in configura- 
tion or energy space, or as characterized by their evolv- 
ing spectral distributions, it is evident that these wildly 
chaotic orbits can, and often do, diffuse more or less freely 
through phase space. 

Unlike the distinction between regularity and chaos, 
the distinction between sticky and wildly chaotic is not 
absolute. If an initially sticky chaotic orbit segment is 
evolved for a sufficiently long period, it will eventually 
transform itself into a wildly chaotic orbit. As discussed 
elsewhere this suggests that, in the presence of the 
driving, the chaotic phase space regions are partitioned 
by 'entropy barriers', corresponding, e.g., to small holes 
through which chaotic orbits eventually effuse. 

Examples of these different types of orbits are exhib- 
ited in Fig. 1. Here the top row corresponds to a regular 
orbit, the middle two rows to sticky chaotic segments, 
and the bottom two to wildly chaotic segments. One 
of these wildly chaotic segments remains bound energet- 
ically {E < 0) for the duration of the integration. The 
other becomes unbound energetically before t = 2048. In 
each case, the leftmost column exhibits the trajectory of 
the orbit as viewed in terms of the phase space coordi- 
nates X and Vx, each recorded at intervals St = 1.0. The 
second column exhibits the mean drift in energy, SE/E, 
as a function of time, and the third exhibits the Fourier 
transform of the resulting time series. The fourth column 
exhibits estimates of the largest Lyapunov exponent 
computed for times < t, and the last column exhibits the 
power spectrum |2;(r2)|. 

As noted already, the type of orbit correlates not only 
with the overall shape of the Fourier spectrum, but with 
the value of the peak frequencies for which quantities like 
|a;(n)| or E{u!) have the most power. This is illustrated 
in Fig. 2, which exhibits scatter plots of X: the largest 
finite time Lyapunov exponent, and n^^\ the frequency 
for which |a:{ri)| has the most power. Each panel was 
generated from the same set of 800 initial conditions. The 
four columns from left to right correspond to increasing 
amplitudes mg = 0, mg = 0.05, mo — 0.1, and mo = 
0.2. The driving frequency luq — 1.3, and the orbital 



structure is dominated by a 2 : 1 resonance. In particular, 
for sticky orbits |a;(ri)| peaks at fi*^^^ — uJo/2 — 0.65. 
(By contrast, as is evident from Fig. 1, |i?(ri)| peaks at 
^ LUo = 1-3.) For aU the regular orbits, fl'^^^ is 
larger than luo/2, the value of which is indicated by the 
vertical dashed line. Alternatively, orbit segments with 
fi(i) < liJq/2 correspond to wildly chaotic orbits. 

For lower amplitude driving, the relative measure of 
chaotic ~ as opposed to regular - orbits and the typi- 
cal size of the largest finite time Lyapunov exponents for 
the chaotic orbits both tend to be smaller. In particular, 
for energies and potentials which admit no chaos in the 
absence of periodic driving, there appears to exist a min- 
imum threshold amplitude, tuq ~ 0.05, for the onset of 
appreciable amounts of chaos. For lower amplitudes, the 
relative measure of sticky chaotic segments - as opposed 
to wildly chaotic segments - also tends to be higher and, 
as will be discussed later, in some cases the time required 
for transitions from sticky to wildly chaotic tends to be 
longer. 

Of especial interest is the case of potentials which, in 
the absence of periodic driving, admit large measures of 
both regular and chaotic orbits. As illustrated in Fig. 3, 
in this case weak driving tends to make some, albeit typ- 
ically not all, the otherwise regular orbits chaotic; and 
the already chaotic orbits tend to remain chaotic. Few 
if any otherwise chaotic orbits are converted to regular. 
However, it is not true that the driving necessarily in- 
creases the values of the largest finite time Lyapunov ex- 
ponents for all the otherwise chaotic orbits. Indeed, for 
sufficiently small mo the driving is comparably likely to 
increase or decrease the value of x- Significantly, though, 
as mo is increased the degree of chaos always appears to 
increase, at least initially. For sufficiently large values of 
mo all the otherwise regular orbits become chaotic and, 
for the otherwise chaotic orbits, there is a systematic in- 
crease in the values of the largest x- However, there will 
still be a few otherwise chaotic orbits which become less 
chaotic in the presence of driving. One other obvious 
trend evident from Fig. 3 is the systematic decrease in 
the values of x(mo^O) at late times. This decrease re- 
flects the fact that the vast majority of the chaotic orbits 
have become wildly chaotic, and that many of these have 
drifted to less negative energies, where the orbital time 
tjj, and hence the Lyapunov time 1/x, is very long. 

Whether, in the presence of periodic driving, an initial 
condition results initially in sticky as opposed to wildly 
chaotic behavior does not appear to correlate strongly 
with its behavior in the absence of driving. For example, 
there is no obvious preferential tendency for otherwise 
regular initial conditions to yield sticky behavior. For 
a fixed potential, the likelihood that an initial condition 
yields a sticky, as opposed to wildly chaotic, orbit is de- 
termined more by the amplitude mo than by the regular 
or chaotic behavior exhibited by the initial condition in 
the absence of periodic driving. When mo is small, sticky 
orbits tend to be more common; as mo increases, the rel- 
ative measure of wildly chaotic segments increases. 
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FIG. 1: Top row: the phase space trajectory of a typical regular orbit evolved in a triaxial Plummer potential with axis 
ratios a?' : \? : (? = 1.25 : 1 : 0.75, with initial energy E = —0.8, mo = 0.1 and wo = 1-3, along with the relative energy 
oscillations, 5E{t)/E, the power spectrum the finite time Lyapunov exponents x computed for times T <t, and the 

power spectrum, |a;(f2)|. Second row: the same for a 'nearly regular' sticky chaotic orbit, with the peak frequency for 
satisfying n'"*^' = ljo and the peak frequency for equally uJo/2. Third row: a more irregular sticky orbit which still has 

the same peak frequencies. Fourth row: a wildly chaotic orbit with spectra peaked at lower frequencies. Fifth row: a more 
wildly chciotic orbit that rapidly diffuses toward larger radii and lower frequencies. 
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FIG. 2: Finite time Lyapunov exponents as a function of 
peak frequency n*-^' for regular (n'^' > aJo/2), sticky (f2^^^ w 
aJo/2), and wildly chaotic (fi*-^^ < loo/2) orbits, generated 
from a microcanonical sampling of 800 initial conditions with 
energy E = —0.8, evolved in a triaxial Plummcr potential 
with 0? : : (? = 1.25 : 1 : 0.75 and driving frequency 
a; = 1.3, integrated for times t = 1024, 2048, 4096, 8192, and 
16384, with amplitudes: (a)-(f) mo = 0, (g)-(l) mo = 0.05, 
(m)-(r) mo = 0.1, (s)-(x) mo = 0.2. 
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FIG. 3: Finite time Lyapunov exponents for driven versus 
unperturbed orbits, integrated from the same microcanoni- 
cal set of 800 initial conditions with E = —0.8 in a triaxial 
Dehnen potential with axis ratios : 6^ : = 1.25 : 1 : 0.75 
and driving frequency u) = 2.75 at varying integration times 
t = 1024, 2048, 4096, 8192, and 16384 for: (a) - (e) mo = 0.1, 
(f) - (j) mo = 0.2 and (k) - (o) mo = 0.3. 



The relative abundance of different orbit types does 
not appear to exhibit a simple, universal dependence on 
shape. In some cases, there is more chaos for parameter 
values where the equipotcntial surfaces arc oblate and/or 
prolate, rather than spherical; in others spherical poten- 
tials exhibit more chaos. The one trend that does emerge, 
at least for the models considered here, is that genuinely 
triaxial potentials with a ^ c tend to exhibit more 
chaos than potentials with either spherical or spheroidal 
symmetry. 

To quantify at least partially the transition rate be- 
tween sticky and wildly chaotic behavior, one can, for 

an orbit ensemble of fixed initial energy, begin by iden- 
tifying the initial conditions that correspond to chaotic 
orbits, then identifying the range of energies which, for 
these initial conditions, correspond to sticky behavior, 
and finally determining for each orbit generated from a 



chaotic initial condition the first time that it acquires an 

energy outside this range. The obvious quantity of in- 
terest is the fraction N(t)/N{0) of those chaotic orbits 
which have not yet 'escaped' within time t. When, as 
for the example presented below, the orbital structure is 
dominated by a 2 : 1 resonance, regular orbits can be re- 
jected immediately as having 0^^' > ujq/2. However, the 
computed fraction does depend somewhat on the precise 
choice of energy range, which is not completely obvious. 
Fortunately, though, as discussed below, modest changes 
do not significantly alter the final results. 

Examples of such an analysis, all generated for the 
same ensemble, are exhibited in Fig. 4, each panel of 
which plots N{t)/N{0), the relative number of chaotic 
orbits which are sticky, as a function of time t for differ- 
ent choices of amplitude toq. The top two panels are for 
oblate (a^ = 62 = 1^^2^ o.75) and triaxial (a^ = 1.25, 
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FIG. 4: (a) Fraction N{t)/N{0) of chaotic orbits that remain 
'sticky' after time t, computed for ensembles of 1200 initial 
conditions with energy E — —0.8, evolved in an oblate Plum- 
mer potential with axis ratio : : — 1 : 1 : 0.75, pulsed 
with frequency loq = 1.3 and variable amplitude mo. Solid 
curves with increasing thickness correspond to increasing am- 
plitudes mo = 0.15, 0.2, 0.25, 0.3. Dashed lines represent best 
linear fits, (b) The same for a triaxial Plummer potential 
with : 6^ : = 1.25 : 1 : 0.75, again with E = -0.8 and 
ujo ~ 1.3. (c) The same for an oblate Dehnen potential with 

: 6^ : = 1 : 1 : 0.75, again with E — —0.8 but now with 
uJo = 2.75. (d) The same for a triaxial Dehnen potential with 

: 6^ : = 1.25 : 1 : 0.75 with E = -0.8 and cuo = 2.75. 

— 1, (? — 0.75) Plummer potentials which, in the ab- 
sence of driving, admit few if any chaotic orbits. The 
lower panels are for oblate and triaxial Dehnen poten- 
tials with the same axis ratios, which admit a consid- 
erable amount of chaos even in the absence of driving. 
From this Figure two points are immediately evident. 
The first is that, in every case, after a very rapid initial 
decrease, In N decreases systematically in a fashion that 
is roughly linear in time. This implies that, at least in 
a first approximation, transitions from sticky to wildly 
chaotic can be viewed a Poisson process, 

7V(t) =iVoexp(-Ai), (5) 

with A^o the number of chaotic orbits which are initially 
sticky and A the escape rate. This exponential decrease 
is the same behavior observed for diffusion through can- 
tori in two-degree-of-freedom systems (see, e.g., |9|), thus 
reinforcing again the notion that these transitions involve 
passage through an entropy barrier. 

The other striking point is that the dependence of A on 
mo is very different for the Plummer and Dehnen models. 
For the Plummer models A is rather nearly independent 
of driving amplitude; but, for the Dehnen models A in- 
creases significantly as mg increases. This is an example 



of a more general result: If, in the absence of driving, 
the orbits are all regular, A exhibits little if any depen- 
dence on mg. If, however, a significant fraction of the 
orbits are already chaotic in the absence of driving, A 
is a strongly increasing function of mp. The degree of 
chaos exhibited by a model, as probed by the relative 
measure of chaotic orbits or the size of a typical finite 
time Lyapunov exponent, need not correlate especially 
with the amount of chaos that is present in the absence 
of driving. However, the presence or absence of chaos in 
the unperturbed models does impact the transition rate 
between sticky and wildly chaotic behavior. 

But why do plots of N{t) show an initial much more 
rapid decrease, and how sensitive are these results to the 
energy range associated with sticky orbit segments? The 
basic conclusion is that increasing the size of the sticky 
range slows the early time rapid decrease, but has only 
a minimal impact on the computed transition rate A. 
The results exhibited in Fig. 4 were computed assuming 
a range — 1.5mo < E < Eq -\- 1.5mo, with Eq the 
initial unperturbed energy and mg the amplitude of the 
perturbation, but, in most cases, the best fit A changes 
only by 10% or so if 1.5 is increased to 3.0 16]. The 
validity of this criterion was also tested by comparing the 
number of orbits that this prescription classifies as sticky 
at time t = 4096 with the number of orbits which still 
appear locked to the driving frequency with $7^^^ = uiq/2. 

The obvious interpretation is that the initial rapid de- 
crease reflects wildly chaotic orbits which require a finite 
time to escape from the sticky range: if the range is made 
larger, the time required for escape becomes longer. This 
suggests in turn that the most accurate estimate of the 
number of sticky orbits involves assuming the smallest 
range which clearly includes all sticky orbits, the crite- 
rion used to justify a value « 1.5. Given such a choice, 
one also appears justified in interpreting A^o in cq. (2.5) 
as the relative fraction of the chaotic orbits which are 
sticky at or near t — 0. 

Given this interpretation of A'o, one can easily esti- 
mate the relative abundance of regular, sticky, and wildly 
chaotic orbits at time t — 0. The results for the en- 
sembles used to construct Fig. 4 are exhibited in Fig. 5. 
It is evident that, in each case, the relative number of 
chaotic orbits is an increasing function of mo, but the 
relative number of initially sticky orbits can vary signifi- 
cantly. For the Plummer potentials, the relative measure 
of sticky orbits tends to be very small irrespective of the 
magnitude of mo. For the Dehnen potentials, however, 
sticky orbits are comparatively common, especially for 
smaller amplitude driving. In particular, for mo < 0.15 
there are more sticky than wildly chaotic orbits for the 
oblate Dehnen potential. It should, however, be stressed 
that the relative measure of sticky versus wildly chaotic 
model need not increase monotonically as a potential is 
made less symmetric. For example, ensembles evolved in 
a spherical Plummer potential perturbed with the same 
amplitude mg tend to exhibit substantially larger relative 
measures of sticky orbits. 
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FIG. 5: (a) Relative number of chaotic orbits (solid line and 
diamonds), sticky chaotic orbits (dashes and squares), and 
wildly chaotic orbits (dot-dashed lines and triangles) ai t = 0, 
plotted as a function of mo, for the ensembles in Fig. 4, as 
predicted by the best linear fits. 



9V dx' 



(6) 



Here rjdx^/dt represents a dynamical friction and is a 
Gaussian-distributed random force which, for specificity, 
is assumed to sample an Ornstein-Uhlenbeck process. 

In the context of a galaxy or charged-particle beam, 
random perturbations associated with an external en- 
vironment could correspond [Tll | to extrinsic noise, for 
which 77 = 0. Alternatively, effects associated with grain- 
iness or internal substructures could correspond to intrin- 
sic noise, with a nonzero i] related to the moments of 
by the Fluctuation-Dissipation Theorem. In particular, 
the effects of graininess can be modeled (see, e.g., ^3]) as 
intrinsic white noise, which can be viewed as a singular 
limit of an Ornstein-Uhlenbeck process with vanishing 
autocorrelation time. 

The assumption that F^ is Gaussian noise sampling 
an Ornstein-Uhlenbeck process implies that the stochas- 
tic process is characterized completely by its first two 
moments, 

{F\t))=0 and 



Assuming that these conclusions are robust, they have 
obvious implications for chaotic phase mixing in systems 
subjected to periodic perturbations. Relatively low am- 
plitude periodic driving can suffice to trigger large mea- 
sures of chaotic orbits. However, these chaotic orbits 
are apt to be primarily sticky and, as such, unlikely to 
exhibit large-scale secular changes in energy. In the con- 
text of charged-particle beams, this would suggest that 
periodic driving could facilitate a comparatively rapid 
loss of coherence, but that it would not necessarily trig- 
ger the formation of a halo. If, however, the amplitude 
of the driving is increased, a larger measure of the or- 
bits will be wildly chaotic; and a substantial fraction 
of these could drift to higher energy states associated 
with an extended halo. This expectation is consistent 
with real experiments involving charged-particle beams 
(e.g., 17|), which indicate that particle bunches which 
undergo larger amplitude oscillations (larger initial 'mis- 
match') tend to form more extended halos. 



III. LOW AMPLITUDE NOISE 

The objective now is to describe how distinctions 
amongst regular, sticky chaotic, and wildly chaotic be- 
havior and, especially, the rate of transitions between 
sticky and wildly chaotic behavior, are impacted by the 
introduction of weak perturbations, modeled as friction 
and/or Gaussian noise. The conclusions derive from an 
analysis of orbits generated as solutions to the perturbed 



{F^{t^)F^m=S^' K{\t,-t2\), (7) 

where ( . ) denotes a statistical average and the autocor- 
relation function K satisfies 

X(|ii -t2\)= A' exp(-|ti - t2|Ac). (8) 

Here A represents the typical 'size' of the random force 
and the autocorrelation time tc the time scale over which 
it changes appreciably. The limit tc ^ corresponds to 
white noise. State-independent, additive noise has 77 and 
K independent of phase space coordinates. Multiplica- 
tive noise arises when these quantities are functions of 
and/or . 

Three different variants were considered: 

1. Additive noise, both white and colored, without fric- 
tion. 

2. Additive noise with friction related by the Fluctua- 
tion-Dissipation Theorem, 

/+00 
dTK{T) = 2eri = D, (9) 
-00 

where O represents a temperature and D the diffusion 
constant entering into a Fokker-Plank description (111] . 
(This implies A^ — Qrj/tc.) On physical grounds, O was 
typically selected to be comparable in magnitude to the 
initial energy of the orbit (s) in question. 

3. Multiplicative noise with K{y^) satisfying 

K{v')^Kovy{v''), (10) 
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with Kq the autocorrelation function appropriate for ad- 
ditive noise, v'^ the squared speed of the orbit, and (w^) 
the mean squared speed of an orbit with the specified 
initial energy. 

There is no uniformly accepted definition of chaos or 
Lyapunov exponent for stochastic differential equations. 
However, the following prescription, implemented in this 
paper, appears to capture one's physical intuition in 
terms, e.g., of the visual appearance of an orbit. Specif- 
ically, in the context of this paper x will be defined in 
the usual way by comparing the evolution of two nearby 
initial conditions, assuming, however, that both orbits feel 
exactly the same noise and, for consistency, the same fric- 
tion. This entails solving a perturbation equation 



dx^dx^ 



5x-' , 



(11) 



where the Hessian is evaluated along the original noisy 
orbit. This prescription obviates the necessity of deciding 
how to relate concrete realizations of a stochastic process 
at two nearby phase space points. 

Consider first the effects of white noise. If, in the ab- 
sence of noise, most of the orbits are not already chaotic, 
the most conspicuous effect resulting from the addition of 
noise is to increase the relative abundance of chaotic or- 
bits. This is, for example, illustrated in Fig. 6, generated 
for the same Plummer models used to produce Fig. 4, 
which exhibits the relative measure fc of chaotic orbits 
as a function of amplitude ry. For all three choices of mo, 
the noise has a discernible effect even for amplitudes as 
small as 77 = 10^''; for 77 as large as 10~^ almost all the 
orbits have become chaotic. Interestingly, however, the 
noise does not have a significant impact on the size of 
a typical finite time Lyapunov exponent for those orbits 
which are chaotic. 

The empirical distinction between sticky and wildly 
chaotic orbits persists in the presence of relatively low 
amplitude noise, and transitions between sticky and 
wildly chaotic still appear to sample a (near-)Poisson pro- 
cess. However, as was the case in the absence of noise, 
the observed behavior of orbit ensembles depends consid- 
erably on whether or not the orbits are almost all chaotic. 

For systems like the Plummer potentials, where there 
is little chaos in the absence of driving, the rate A associ- 
ated with transitions between sticky and wildly chaotic is 
largely independent of 77. Increasing the noise amplitude, 
which increases the overall abundance of chaotic orbits, 
also increases the abundance of initially sticky chaotic 
segments; but the rate at which these are transformed 
into wildly chaotic orbits is largely unaffected unless the 
noise becomes very large in amplitude. 

Alternatively, for systems like the Dehnen potentials, 
where there is a good deal of chaos even in the ab- 
sence of driving, the addition of noise definitely accel- 
erates transitions from sticky to wildly chaotic. This is, 
e.g., illustrated in the left hand column of Fig. 7, which 
exhibits N{t)/N{0), the relative number of orbit that 
are still sticky, for the triaxial Dehnen potential used to 
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FIG. 6; Relative measure of chaotic orbits, fc, computed 
for a microcanonical ensemble of 1200 initial conditions with 
energy E — —0.8 in a pulsating triaxial Plummer potential 



with axis ratio a 



1.25 



0.75, evolved with 



ujo ~ 1-3, mo — 0.1 (empty circles with solid line), mo — 0.15 
(crosses with short-dashed line), mo — 0.2 (asterisks with 
long-dashed line) and mo = 0.3 (filled circles with dash-dotted 
line) in the presence of additive white noise with Q = 1.0 and 
variable amplitude 77. 



generate Fig. 4. The qualitative behavior manifested in 
these panels is very similar to that associated with dif- 
fusion through cantori in two-degree-of-freedom Hamil- 
tonian systems (c/. 9;]), where one can argue that, by 
wiggling orbits, 'noise helps orbits find holes' in entropy 
barriers. 

The situation for larger amplitude noise is rather differ- 
ent. Seemingly independent of the form of the potential, 
for 77 > 10^'^ — 10~^ escapes happen so quickly that the 
distinction between sticky and wildly chaotic becomes 
essentially meaningless. This can be understood by ob- 
serving that, when the noise has such a large amplitude, 
it will pump energy into the orbit so rapidly that, even 
if the driving has a relatively minimal effect, one cannot 
think of the orbit as being restricted to a nearly constant 
energy hypersurface. 

Although these conclusions were derived allowing only 
for additive white noise in the absence of any dynamical 
friction, they appear to be comparatively insensitive to 
the form of the perturbation. For example, the presence 
or absence of a dynamical friction of comparable ampli- 
tude 77 is largely irrelevant. Moreover, there is no indica- 
tion that allowing for 'generic' multiplicative noise would 
alter the basic picture. For example, allowing for mul- 
tiplicative noise with an autocorrelation function scaling 
as in eq. (3.5) has only a minimal effect on the relative 
number of chaotic orbits or the transition rate between 
sticky and wildly chaotic behavior. All that really seems 
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FIG. 7: (a) Fraction N(t)/N{0) of chaotic orbits that remain 
'sticky' after time t, computed for ensembles of 1200 initial 
conditions with E = —0.8, evolved in a triaxial Dehnen po- 
tential with axis ratio a? : : c? = 1.25 : 1 : 0.75, ujq = 2.75, 
and mo = 0.1, allowing for additive white noise with variable 
amplitude. Dashed lines correspond to integrations without 
white noise, while the magnitude of 77 increases with increas- 
ing line thickness through values rj — 10~^, 10~®, 10~^, 10~^. 
In each case © = 1. (b) The mean energy {E{t)) computed 
for the same ensemble, (c) and (d) The same for mo — 0.15. 
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FIG. 8: Mean value of the quantity AV^, as defined in ll.SI 
computed for 1200 orbits with initial energy E = —0.8 evolved 
in a pulsating triaxial Plummer potential with axis ratio : 
: (? — 1.25 : 1 : 0.75, wq = 1.3, and variable amplitude 
(a) mo = 0.1 and (b) mo = 0.2. Dashed lines correspond 
to integrations without white noise, while the magnitude of 
white noise r) increases with increasing line thickness through 
values r) = 10-^ 10"®, 10-^ lO"*. (c) - (d) The same for a 
triaxial Dehnen potential with axis ratio 1.25 : 1 : 0.75, and 
uo = 2.75. 



to matter is the amplitude of the noise. Moreover, even 
this dependence on amplitude is comparatively weak. In- 
spection of Figs. 6 and 7 indicates that there is only a 
roughly logarithmic dependence on rj. 

Also of interest is how the combined effects of periodic 
driving and friction/noise impact the mean energies of 
typical particle orbits. Examples of such behavior, again 
computed for orbits in a pulsed, triaxial Dehnen potential 
perturbed only by additive white noise, are exhibited in 
the right hand column of Fig. 7. 

But how are such data to be interpreted? It follows 
immediately from the Langevin equation (3.1) that, al- 
lowing only for periodic driving and additive white noise, 
the mean energy {E) associated with some ensemble will 
satisfy 

d(E) /dV\ 

where {dV/ dt) represents an average computed along the 
orbits. The noise leads directly to a linear increase in 
energy at a constant rate 30ry, independent of the form of 
the driving; the driving occasions an additional increase 
which will in general depend on the form of the perturbed 
orbits. Given the mean value {E{t)) computed for orbits 
that are subjected to both noise and periodic driving, it 
is thus natural to subtract off the linear increase 30r/i 
directly attributable and to compare the residual with 



the mean {E{t)) computed in the absence of noise. To 
the extent that the former quantity is larger, one can 
infer that noise enhances the energy diffusion associated 
with periodic driving in the sense that the mean change 
in energy is larger than what would obtain allowing only 
for the separate effects of driving and noise. 

The results of such an analysis are exhibited in Fig. 8, 
each panel of which plots the quantity 

AV = J (dV/dt) dt, (13) 

derived for the same set of initial conditions, evolved in 
the same potential with the same amplitude toq and same 
luq but in the presence of additive white noise with dif- 
ferent amplitudes rj. For both the Plummer and Dehnen 
potentials, it is evident that, at least for low amplitudes 
niQ, increasing 77 tends to increase AV, i.e., that the noise 
and driving reinforce one another in a superadditive fash- 
ion. However, for higher amplitudes mo, more variability 
is observed, the one uniform feature being that the rela- 
tive degree of reinforcement is less. In this sense, one can 
say that noise is more important as a source of energy 
diffusion in systems subjected to relatively low amplitude 
driving. 

But how is all this effected by allowing for colored noise 
with a finite autocorrelation time? If the autocorrelation 
time tc is sufficiently short compared with the dynamical 
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time, colored noise has virtually the same effect as white 
noise. However, when tc becomes comparable to the or- 
bital time to the effects of the noise become weaker; and, 
for tc sufficiently large compared with to the noise has 
only a relatively minimal effect. 

Figure 9 illustrates how, for the same initial conditions 
used to generate Fig. 6, the relative number /c of chaotic 
orbits depends on both amplitude r] and autocorrelation 
time tc- The top panel illustrates clearly that, for fixed 
r^, increasing reduces the degree to which the noise 
converts regular orbits to chaotic. However, it is also ap- 
parent that, for sufficiently large 77, almost all the orbits 
still become chaotic. The autocorrelation time must as- 
sume a value even larger than tc — 32 if the effects of the 
noise are to 'turn off' almost completely for amplitudes 
as large as 77 = 10~^. 

The panel a) of Fig. 10, generated for orbits in tri- 
axial Plummer, demonstrate that the rates of transition 
between sticky and wildly chaotic behavior are largely 
independent of the autocorrelation time tc- The panel 
c) of Fig. 10, representing orbits in triaxial Dehnen po- 
tentials, show that increasing the autocorrelation time tc 
tends also to weaken the effects of noise in accelerating 
transitions between sticky and wildly chaotic behavior. 

The panels b) and d) of Fig. 10, generated for the same 
initial conditions, demonstrates that increasing autocor- 
relation time tc decreases the rates of the mean change 
in energy. For both potentials, the curve for tc = 0.5 
coincides very nearly with the white noise curve, and the 
appreciable decreases in {E{t)) arise for tc > 1.0. 



IV. FREQUENCY VARIATIONS 

The objective now is to explore how the results de- 
scribed in the preceding Section are altered if one al- 
lows instead for variations in the pulsation frequency, 
again modeled as Gaussian noise sampling an Ornstein- 
Uhlenbeck process. What this entails is allowing for a 
variable frequency u){t) = cjq + 6u){t), where the random 
frequency shift Su) is characterized by moments 



and 



{Soj{t)) = 



{d(j{ti)d(j{t2)) = exp(-|ti - t2\/tc 



(14) 



This implies in particular that = + "^he 

form of the experiments was the same as in the preced- 
ing Section. The computations were performed assuming 
F'' = rj = 0; i.e no ordinary noise. 

As discussed in the preceding Section, allowing for 'or- 
dinary' white noise tends to increase the relative abun- 
dance of chaotic orbits, especially for systems which, in 
the absence of driving, admit few if any chaotic orbits. 
Allowing for a 'noisy' frequency with vanishing autocor- 
relation time has exactly the opposite effect. For sys- 
tems where, in the absence of driving, little if any chaos 
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FIG. 9: Relative measure of ehaotie orbits, fc, computed 
for a microcanonieal ensemble of 800 initial conditions with 
energy E = —0.8 in a pulsating triaxial Plummer potential 
with axis ratio c? ■.}?■.(? = 1.25 : 1 : 0.75, evolved with 
frequency ojo = 1-3 and mo = 0.1 in the presence of noise 
with varying amplitude r\ and autocorrelation time ta- In (a): 
solid line with empty circles represents white noise (limit tc 
0), short-dashed line with crosses tc ~ 0.5, long-dashed line 
with squares tc — 2, and dash-dot line with triangles tc = 8, 
dash-dot-dot-dot line with filled circles tc = 32. In (b): solid 
line with empty circles represents rj = 10~^, short-dashed 
line with crosses rj = 10~®, long-dashed line with squares 
ri = 10~®, dash-dot line with triangles r; — 10"'', dash-dot- 
dot-dot line with filled circles r/ = 10^'*. Horizontal dashed 
line corresponds to the fraction of chciotic orbits for the same 
ensemble evolved in the absence of any noise. 
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FIG. 10; (a) Fraction N{t) /N{0) of chaotic orbits that remain 
'sticky' after time t, computed for ensembles of 800 initial 
conditions with E = —0.8, evolved in a triaxial Plummer 
potential with axis ratio : : = 1.25 : 1 : 0.75, wo = 1.3, 
and mo = 0.1, allowing for colored noise with r] = 10~* and 
variable tc- The magnitude of tc increases with increasing 
line thickness through values tc — 0.5, 1, 2, and 4. The long- 
dashed line corresponds to white noise. The short-dashed line 
corresponds to the same ensemble evolved in the absence of 
any kind of noise, (b) The mean energy {E{t)) computed 
for the same ensembles, (c) and (d) The same for ensembles 
evolved in a triaxial Dehnen potential with : : = 1.25 : 
1 : 0.75 and ojo = 2.75, again with mo = 0.1 and rj = 10~^. 



is present, making the frequency 'stutter' by introduc- 
ing even comparatively low amplitude white noise im- 
pulses tends instead to decrease the relative abundance 
of chaotic orbits. This is, e.g., evident from Fig. 11, 
which was generated for the same initial conditions as 
Fig. 6, now allowing for white noise variations in u) with 
different amplitudes 0^. 

As is illustrated in Fig. 12, allowing for random varia- 
tions in frequency tends also to suppress transitions from 
sticky to wildly chaotic behavior and to decrease the 
amount of energy that the driving can pump into the 
orbits. White noise variations in frequency dramatically 
reduce the effects of driving as a source of energy diffu- 
sion. 

These differences are reflected by the fact that even 
very weak perturbations, with amplitudes as small as 
Q.^ = 10~^, completely alter the Fourier spectra of the 
orbits. Examples thereof are provided in Fig. 13 for both 
triaxial Plummer and triaxial Dehnen potentials. In each 
case, it is evident that, in the absence of the perturba- 
tion, there is a dominant resonance, = Uo/2 = 0.65 
for the Plummer potential and O^^^ = wo/2 = 1.375 for 
the Dehnen potential. But, in each case, it is also appar- 
ent that, even for noise as weak as = 10^^, the effects 
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FIG. 11: Relative measure of chaotic orbits, /c, computed 
for a microcanonical ensemble of 800 initial conditions with 
energy E = —0.8 in a pulsating triaxial Plummer potential 
with axis ratio a'' ■.}?■.(? = 1.25 : 1 : 0.75, evolved with 
oJo = 1-3, mo = 0.1 (empty circles with solid line), mo = 0.15 
(crosses with short-dashed line), mo = 0.2 (asterisks with 
long-dashed line) and mo = 0.3 (filled circles with dash-dotted 
line), allowing for 'white noise' random variations in frequency 
with variable amplitude f2^. 



of this resonance have entirely disappeared. 

Also evident from this Figure is the fact that even a 
very small suffices to prevent orbits from diffusing to 
very large energy. In the absence of frequency variations, 
both potentials admit orbits that diffuse towards very 
large radii and, consequently, have very small peak fre- 
quencies. However, for ^ > 10^"'', there arc no orbits 
with fi^^^ < 0.4. This does not, however, necessarily im- 
ply that there are no chaotic orbits at all. For the case of 
the Dehnen potential there arc still significant numbers 
of chaotic orbits with larger values of Q^^^ 

One might naively have expected that the stuttering 
simply 'turns off' the effects of the driving, so that or- 
bits pulsed with mo = 0.1 and = 10"^, lO'^, or 10"^ 
behave very much like orbits evolved with mo = fi^ = 0. 
This, however, is not the case. It is, for example, clear 
that, for the Dehnen model, the distribution of peak 
frequencies for an unpulsed, constant frequency integra- 
tion has far more structure than the distributions for 
the pulsed, variable frequency models. Even though fre- 
quency variations can dramatically reduce the efficacy of 
the driving as a source of energy diffusion, they do not 
transform orbits back to the form they had in the absence 
of driving. 

Adding color strengthens the effects of the perturba- 
tion. An example of this tendency is illustrated in Fig. 14, 

which exhibits the fraction of chaotic orbits for the en- 
sembles of 800 orbits in a pulsed Plummer potential, al- 
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FIG. 12: Fraction N{t)/N{0) of chaotic orbits that remain 
'sticky' after time t, computed for ensembles of 800 initial 
conditions with E = —0.8, evolved in a triaxial Dehnen po- 
tential with axis ratio a'' :}?■.(? = 1.25 : 1 : 0.75, wo = 2.75, 
and mo = 0.1, allowing for 'white noise' variations in fre- 
quency with variable amplitude n. The magnitude of the 
noise increases with increasing line thickness through values 
= 10-^ 10-*^, 10-^ and 10"*. The dashed hric corre- 
sponds to integrations without frequency variations, (b) The 
mean energy {E{t)) computed for the same ensembles, (c) 
and (d) The same for mo = 0.2. 



lowing for random variations in frequency with different 
values of and t,.- It is clear from the bottom panel 
that increasing the value of autocorrelation time from the 
tc ^ white noise limit results in a systematic increase 
in the relative measure of chaotic orbits. Also evident 
is the fact that even a relatively small value tc = 0.5 is 
sufficient to significantly increase the number of chaotic 
orbits relative to the white noise limit. As is evident from 
the lower panel, in the long autocorrelation time limit, 
— > oo, the effects of the colored noise gradually turn 
off. This is to be expected, if one views tc ^ oo as an 
adiabatic limit. The time-scale on which this turning off 
occurs is directly proportional to the amplitude of col- 
ored noise, i.e, it takes longer autocorrelation times for 
the effects of the stronger noise to reach the adiabatic 
limit. 

The left panel of Fig. 15 illustrates the fact that allow- 
ing for a finite autocorrelation time enhances the transi- 
tions from sticky to wildly chaotic behavior. The right 
hand panel illustrates the fact that, for fixed amplitude 
fl^, integrations with nonvanishing autocorrelation time 
tc will pump more energy into the orbits than correspond- 
ing integrations with white noise or in the absence of 
noise. The observed approach of the thick solid lines 
representing large autocorrelation time tc to the short- 
dashed lines corresponding to integrations in the absence 



of any noise illustrates that the effects of the noise are 
turning off in the tc oo limit. 

One other point would seem significant. For the case 
of 'ordinary' noise, incorporated as a random force in 
the equations of motion, typically a comparatively large 
autocorrelation time tc ~ 1 is required before one sees 
appreciable differences from the white noise limit. By 
contrast, however, for the case of 'noisy' variations in fre- 
quency even an autocorrelation time as small as tc = 0.5 
yields substantial differences from the white noise limit. 



V. CONCLUSIONS 

Orbits in three-dimensional potentials subjected to pe- 
riodic driving, V{x^,t) = (1 -|- mo sina;t)Vb(a;*), divide 
into two distinct types, between which transitions are 
seemingly impossible: regular orbits, which are strictly 
periodic and have vanishing Lyapunov exponents, and 
chaotic orbits, which are aperiodic and have at least one 
positive Lyapunov exponent. Viewed over finite time in- 
tervals, the chaotic orbits divide in turn into two types: 
'sticky' chaotic segments, which are locked to the driving 
frequency and exhibit little systematic energy diffusion, 
and 'wildly chaotic' segments, which are not so locked 
and can exhibit significant energy diffusion. The latter 
distinction is not absolute, transitions between sticky and 
wildly chaotic apparently involving motion through en- 
tropy barriers. 

The relative abundance of sticky versus wildly chaotic 
segments depends considerably on the form of the po- 
tential. However, sticky segments tend overall to be 
more common for low amplitude driving, i.e., the rel- 
ative phase space measure associated with sticky orbits 
tends to be larger. Increasing the amplitude also tends 
to increase the transition rate between sticky and wildly 
chaotic behavior, which suggests that the 'holes' associ- 
ated with the entropy barriers become larger as mo in- 
creases. 

The properties of orbits in such time-periodic poten- 
tials can be impacted significantly even by very low am- 
plitude perturbations, both random impulses added to 
the equations of motion and random variations in the 
driving frequency, in each case modeled as Gaussian noise 
sampling an Ornstein-Uhlenbeck process. For fixed am- 
plitude of the perturbation, the largest effects arise in the 
white noise limit of vanishing autocorrelation time tc and 
decrease systematically with increasing tc- The details of 
the perturbation tend to be comparatively unimportant. 
For example, for the case of 'ordinary' noise the presence 
or absence of friction of comparalile magnitude and/or 
making the noise multiplicative as opposed to additive 
has a comparatively minimal effect. Even the precise 
strength of the perturbation is relatively unimportant, 
the response exhibiting only a weak, roughly logarithmic 
dependence on amplitude. 

Ordinary noise tends to enhance the effects of periodic 
driving by increasing the relative measure of chaotic or- 
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bits, especially in potentials which contain little if any 
chaos in the absence of driving. It also acts to acceler- 
ate transitions between sticky and wildly chaotic behav- 
ior, thus facilitating more efficient energy diffusion. The 
effects of noise in frequencies of pulsations have more 
complex effects. White noise suppresses the effects of 
periodic driving by decreasing transitions between sticky 
and wildly chaotic behavior. Colored noise fluctuations, 
at least for autocorrelation times on the order of the dy- 
namical time, enhances the effects of periodic driving. 
These effects gradually turn off in the adiabatic t,, ^ oo 
limit. 

Ordinary noise appears to act by enhancing the effects 
of the resonances associated with the periodic driving. 
Wiggling orbital trajectories can allow an otherwise reg- 
ular orbit to be impacted by a resonance which it might 
not otherwise 'feel', thus causing it to become chaotic. 
And similarly, the accelerated rate of transitions between 
sticky and wildly chaotic behavior can be attributed to 
the fact that such wiggling helps already chaotic orbits 
'find' holes in the entropy barrier. Indeed, the effect of 
noise in enhancing such transitions is very similar qualita- 
tively to the effect of noise as a source of accelerated dif- 



fusion through cantori in time-independent two-degree- 
of- freedom Hamiltonian systems @. 

Random white noise variations in the driving frequency 
appear to weaken the resonances. Making the frequency 
'stutter' by adding even a very low amplitude, high fre- 
quency component to the driving significantly alters the 
Fourier spectra of orbits, suppressing structure which, 
in the absence of the perturbation, indicates the impor- 
tance of the resonances. Variations characterized by non- 
vanishing autocorrelation time enhances the resonances 
because the Fourier spectrum of the frequency now has 
another strong component with which the fundamental 
frequencies of orbits can resonate. 
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APPENDIX A: IMPLEMENTING WHITE AND 
COLORED NOISE 



The second term is transient, so that for t^tc. 



Here we outhne the implementation of noisy equations 
of motion © . The term {t) is a Gaussian-distributed 
random force which samples Ornstein-Uhlenbeck process. 
It has a autocorrelation time tc, which in the limit tc — > 
samples ordinary white noise. F^{t) is a normalized 
solution to the free particle Langevin equation 



du^ 1 , 1 



(Al) 



where {t) is Gaussian-distributed white noise with zero 
mean and unit variance. It is well known that such a vari- 
able is easily generated as a sum of 12 random numbers 
in the range [—1/2,1/2]. The solution to (|A1|) is given 

by m 

u\t) = ui exp{-t/tc) 

t 

+ ^ exp(-i/ic) / exp(C/tc) d^{^2) 



which, after some straightforward algebra, yields 

i^'^) = 7^ + f^o*' - ^) expi-2t/Q. (A3) 



2tc 



2tc 



1 

2L. 



(A4) 



This means that the random variable y/2t^u^ has auto- 
correlation time tc, zero mean and unit variance. 

We determine the proper amplitude of colored noise A 
after substituting F^{t) — y/Au^lt), into the equation ^ 
and equating it with the equation (jSJ: 



{Fm,)F^{t2)) = 



— exp(-|ti - iai/ic) 

A2exp(-|ti-t2|/ic) 
B 



= --exp(-|ti-i2|/ic), (A5) 



from which it is clear that A = D = 2Qri. There- 
fore, the properly scaled noise noise factors are F'^{t) = 
\/2Qrju^{t), where u^{t) is the solution of lAip for colored 
noise, and F'^{t) — ^2Qrj w^{t), where w'^{t) is Gaussian- 
distributed white noise with zero mean and unit variance. 
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FIG. 13: Finite time Lyapunov exponent as a function of 
frequency fi-'^' for 800 orbits with initial energy E = —0.8 in 
a triaxial Plummer potential with axis ratio : : = 
1.25 : 1 : 0.75, pulsating with amplitude mo = 0.1 and loo = 
1.3, allowing for 'white noise' variations with (a) f2 = 0, (b) 
= 10"'^, (c) n = 10"^ and (d) ^ = 10"^ (e) The same 
ensemble evolved in the absence of both noise and driving, (f ) 
- (j) The same for orbits in a Dehnen potential with uiq = 2.75, 
again with : 6^ : = 1.25 : 1 : 0.75 and mo = 0.1. 
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FIG. 14; Relative measure of chaotic orbits, fc, computed 
for a microcanonical ensemble of 800 initial conditions with 
energy E = —0.8 in a pulsating triaxial Plummer potential 
with axis ratio c? ■.}?■.<? = 1.25 : 1 : 0.75, evolved with 
frequency wo = 1.3 and mo = 0.1, allowing for 'colored noise' 
variations in frequency with variable amplitude and au- 
tocorrelation time tc- In (a): solid line with empty circles 
represents white noise (limit tc ^ 0), short-dashed line with 
crosses tc ~ 0.5, long-dashed line with squares tc = 2, and 
dash-dot line with triangles tc = 8, dash-dot-dot-dot line with 
filled circles tc = 32. In (b): solid line with empty circles rep- 
resents fi^ = 10~^ , short-dashed line with crosses ft^ = 10~®, 
long-dashed line with squares = 10~^, dash-dot line with 
triangles = 10~*, dash-dot-dot-dot line with filled circles 
fl^ = 10"'^. Horizontal dashed line corresponds to the frac- 
tion of chaotic orbits for the same ensemble evolved in the 
absence of any noise. 
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FIG. 15: (a) Fraction N{t)/N{0) of chaotic orbits in a 800 
orbit ensemble that remain 'sticky' after time t, computed 
for an ensemble of initial conditions with E = —0.8, evolved 
in a pulsating triaxial Dehnen potential with axis ratio : 
6^ : = 1.25 : 1 : 0.75, wo = 2.75 and mo = 0.1, allowing for 
'colored noise' random variations in frequency with = 10~^ 
and variable tc. The magnitude of tc increases with increasing 
line thickness through values tc ~ 0.5, 2, 8, and 32. The long- 
dashed line corresponds to white noise. The short-dashed line 
corresponds to the same ensemble evolved in the absence of 
any kind of noise, (b) The mean energy {E{t)) computed for 
the same ensembles. 



